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CJ . Abstract 

We present a model of the gravitational field based on two sym- 
metric tensors. The equations of motion of test particles are derived: 
Massive particles do not follow a geodesic but massless particles tra- 
jectories are null geodesies of an effective metric. Outside matter, 
the predictions of the model coincide exactly with General Relativity, 
so all classical tests are satisfied. In Cosmology, we get accelerated 
expansion without a cosmological constant. 

General Relativity(GR) works very well at the macroscopic scaIes[T]. Its 
quantization has proved to be difficult, though. It is non renormalizable, 
which prevents its unification with the other forces of Nature. Trying to make 
sense of Quantum OR is the main physical motivation of String Theories 
[2]. Moreover, recent discoveries in Cosmology [3111] has revealed that most 
part of matter is in the form of unknown matter(dark matter, DM) and that 
the dynamics of the expansion of the Universe is governed by a mysterious 
component that accelerates the expansion(dark energy,DE). Although OR is 
able to accommodate both DM and DE, the interpretation of the dark sector 
in terms of fundamental theories of elementary particles is problematic [5]. 
Although some candidates exists that could play the role of DM, none have 
been detected yet. Also, an alternative explanation based on the modification 
of the dynamics for small accelerations cannot be ruled out [6]. 

In GR, DE can be explained if a small cosmological constant (A) is present. 
At the later stages of the evolution of the Universe A will dominate the 
expansion, explaining the acceleration. Such small A is very difficult to 
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generate in Quantum Field Theory (QFT) models, because in this models A 
is the vacuum energy, which is usually very large. 

One of the most important questions in Cosmology and cosmic structure 
formation is to understand the nature of dark energy in the context of a 
fundamental physical theory [T7]. 

In recent years there has been various proposals to explain the observed 
acceleration of the universe. They involve the inclusion of some additional 
field like in quintessence, chameleon, vector dark energy or massive grav- 
ity; Addition of higher order terms in the Einsten-Hilbert action, like f(R) 
theories and Gauss-Bonnet terms; Modification of gravity on large scales by 
introduction of extra dimensions. For a review, see [7j. 

Less widely explored, but interesting possibilities, are the search for non- 
trivial ultraviolet fixed points in gravity (asymptotic safety [9]) and the notion 
of induced gravity[TO]. The first possibility uses exact renormalization-group 
techniques [11] and lattice and numerical techniques such as Lorentzian trian- 
gulation analysis p!2]. Induced gravity proposed that gravitation is a residual 
force produced by other interactions. 

In a recent paper, [13] a two-dimensional field theory model explore the 
emergence of geometry by the spontaneous symmetry breaking of a larger 
symmetry where the metric is absent. Previous work in this direction can be 
found in [H], [15] and [16] . 

In this paper, we wish to present a model of gravitation that is as close as 
possible to classical GR, but could make sense at the quantum level. The 
main observation is that GR is finite on shell at one loop [18]. In [20], [19] 
we presented a type of gauge theories, 6 gauge theories(DGT): The main 
properties of DGT are: 1) The classical equations of motion are satisfied in 
the full Quantum theory 2) They live at one loop. 3) They are obtained 
through the extension of the former symmetry of the model introducing an 
extra symmetry that we call S symmetry, since it is formally obtained as the 
variation of the original symmetry. When we apply this prescription to GR 
we obtain 6 gravity. Quantization of 6 gravity is discussed in [2T| . 

The impact of dark energy on cosmological observations can be expressed 
in terms of a fluid equation of state p = w{R)p, which is to be determined 
studying its influence on the large-scale structure and dynamics of the Uni- 
verse. 

In this paper we follow the same approach. So we will not include the 
matter dynamics, except by demanding that the energy-momentum tensor of 
the matter fluid is covariantly conserved. This is required in order to respect 
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the symmetries of the model. 

The main properties of this model at the classical level are: a)It agrees 
with GR, outside the sources and with adequate boundary conditions. In 
particular, the causal structure of delta gravity in vacuum is the same as in 
General Relativity. So all standard test are satisfied automatically, b) When 
we study the evolution of the Universe, it predicts acceleration without a 
cosmological constant or additional scalar fields. The Universe ends in a Big 
Rip, similar to the scenario considered in [23]. c) The scale factor agrees 
with the standard cosmology at early times and show acceleration only at 
late times. Therefore we expect that density perturbations should not have 
large corrections. 

It should be remarked that S- gravity is not a metric model of gravity 
because massive particles do not move on geodesies. Only massless particles 
move on null geodesies of a linear combination of both tensor fields. 

It was noticed in [20] that the Hamiltonian of delta models is not bounded 
from below. Phantoms cosmological models [22] , [23] also have this property. 
Although it is not clear whether this problem will subsist in a diffeomorphism 
invariant model as delta gravity or not, we mention some ways out of the 
difficulty at the end. 

Definition of Delta gravity In this section we define the action as well 
as the symmetries of the model and derive the equations of motion. 

We use the metric convention of [5]. The action of 6 gravity is: 



the energy-momentum tensor of matter. is the Ricci's tensor and R is 
the curvature scalar of g^^. g^'^ is a two- contravariant tensor under general 
coordinate transformations. 

The action ([1]) is obtained by applying the prescription contained in [201 
[T9] . That is, we add to the action of general relativity, the variation of it and 
consider the variation 5g^y = g^^ as a new field. Similarly, the symmetries we 
write below are obtained as variation of the infinitesimal general coordinate 




+ 



(1) 



Here n 



, ^2 is an arbitrary constant 
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transformations where the variation of the infinitesimal parameter 5^q = 
is the infinitesimal parameter of the new transformation 6. The last term in 
(II]) is needed to implement the condition T.'^'^ = as an equation of motion 
in order to implement the 6 symmetry (|2]) off shell. This term is not needed 
in vacuum. 

Action ([T]) is invariant under the following transformations((5): 

^gpu{x) = ^i^.y + iijj-^ + 9pp^o,i^ + 9iyp^o,i_i + 9pu,p^Q (2) 
^^p. = ~^ip + ^p^o,p + ^pAo 

From now on we will fix the gauge = 0. This gauge preserves general 
coordinate transformations but fixes completely the extra symmetry with 
parameter ^i^j. 

Equations of motion Varying g^^ we get: 

+ ^-{Rr'' - gpur^nn - (v^v.^'^'^)^^ + 

i^^^^ [^99-'^ ,{9p.~9n) = -^^f' (3) 

where S^'' = {V'^^p + W^^^p - U"^Pp).pp, U'^^^p = \ {g^Pi^f^^^ - ^g"'^ 9f,u9'"')] 
Varying gP'^ we get Einstein equation: 

(^Rpu - IgpuR^ + ^T^u = (4) 

Varying we get.T.Pf = 

Covariant derivatives as well as raising and lowering of indices are de- 
fined using g^j^i,. Notice that outside the sources(r^^ = 0), a solution of ([3]) 
is gP"" = Xgp", for a constant A, since gp^ = and i?^^ = 0. We will have 
9'^" = 9^^ 1 assuming that both fields satisfy the same boundary conditions 
far from the sources. But there exists other solutions in the vacuum. A 
simple case is presented in equation (48) of [21]. This solution is interest- 
ing because any finite size body will look point-like if we watch it from far 
away. Studying the motion of massive and massless test particles , using the 
equations we will derive below, we can see that the parameter j3 produces an 
additional gravitational force. We need further studies to ascertain whether 
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this additional gravitational force can be used to understand dark matter or 
not. 

The equation for g'^'^ is of second order in the derivatives. 

Particle motion in the gravitational field We are aware of the pres- 
ence of the gravitational field through its effects on test particles. For this 
reason, here we discuss the coupling of a test particle to a background grav- 
itational field, such that the action of the particle is invariant under 

In 6 gravity we postulate the following action for a test particle: 

S, = -mj dt^-g^^xf^x^ + 4 j ryV^%, + A^'^" + A'^^^) 
where T^^ is the energy momentum tensor of the test particle: 



2^ J ^-ga^x'^xP 

K'2 = /t2/« is a dimensionless constant. 
That is: 

were g^^ = g^u + X^-u + Xi,-^ . Notice that Sp is invariant under ([2]) and 
t-parametrizations. 

From now on we work in the gauge = 0. 

Since far from the sources, we must have free particles in Minkowski 
space, i.e g^i, ~ rj^y^g^^ ~ rj^^, it follows that we are describing the motion of 
a particle of mass m! = m(l + 

Since in vacuum g^" = g'^'^, the equation of motion for test particles is 
the same as Einstein's. Moreover, the equation of motion is independent of 
the mass of the particle. 

In order to include massless particles, we prefer to use the action [23]: 



2 y (^^"^ " ^ ^3i.u + K2g^.u) x^x + 2v-^g^^±c.r^Si \rn \v gxpX xP) 

This action is invariant under reparametrizations: 

x'{t') = x{t); dt'v'it') = dtv{t); t' = t - e{t) (7) 
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The equation of motion for v is: 



m 

Replacing ([8]) into ([6]), we get back ([5]). 

Let us consider first the massive case. Using ([7]) we can fix the gauge 
V = \. Introducing mdt = dr, we get the action: 

Li = l^J dr (l - {g,^ + 4^^^) x^i" + ^ ^ • ^^J^ ^'^'^ + dxpx'x^)^ (9) 

plus the constraint obtained from the equation of motion for v: 

g^.x^'x^ = -1 (10) 
From Ll the equation of motion for massive particles is derived. We define:0^j, = 

d{x'^x'^Q^j,x'^ gai3 + 2x^^2013) l.^-v- •/3-7 ■ ^J, ■ 

^— —X X Q^uX X 5'/37,a XX „ = (H) 

We will discuss the motion of massive particles elsewhere. 
The action for massless particles is: 

Lo = ^jdt {-v~' {g^, + K'^g^,) x^'x') (12) 

In the gauge = 1, we get: 

Lo = -^j dt {g^,^ + K!2gt,u) x'^x" (13) 

plus the equation of motion for v evaluated at = 1: {g^^ + K^g^u) x^x^ = 
So, the massless particle moves in a null geodesic of g^i, = g^i, + K^g^j^y 
Distances and time intervals In this section, we define the measure- 
ment of time and distances in the model. 

In GR the geodesic equation preserves the proper time of the particle 
along the trajectory. Equation ([11]) satisfies the same property: Along the 
trajectory x^x^g^y is const ant. Therefore we define proper time using the 
original metric g^j^, 



dr = —g^j^ydx^^dx^ = y^—goodx{dx^ = 0) (14) 
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Following [25], we consider the motion of light rays along infinitesimally near 
trajectories and f ll4p to get the three dimensional metric: 

7ii = —idij - ——^) (15) 
000 Soo 

That is, we measure proper time using the metric Qfj^^, but the space geometry 
is determined by both metrics. In this model massive particles do not move 
on geodesies of a four dimensional metric. Only massless particles move on 
a null geodesic of q^^^. So, delta gravity is not a metric theory. 

The Newtonian limit The motion of a non relativistic particle in a weak 
static gravitational field is obtained using ^(^j/ = diag (—1 — 2 Ue, 1 — 2 Ue,l — 2 Ue,l — 2 Ue), 
which solves Einstein equations to first order in e if V^t/ = jKp. 

The solution for g^^ is g^^ = diag (^eU,l+ e (^U -2U^ ,1+ e (^U -2U^ ,1+ e (u -2U 

Solving ([3]), to first order in e we get V^f/ = ^Kp. 

To recover the Minkowsky metric far from the sources, p — 0, we must 
require there:t/ — )■ 0, [7 — )■ — e~^. 

dnD implies ^ = -0,^ with <p = U - n'^{2U + tj). 

The Newtonian potential satisfies = 1(1 — 3^2)/}, 1^2! <^ 1. The 
whole effect is a small redefinition of Newton constant. 

Gravitational red shift experiments can be used to put bounds on k^- 
According to ([H]), the shift in frequency of a source located at Xi, compared 
to the same source located at X2 due to the change in gravitational potential 
is: "'^'^^ = 4>n{.X2) — 4>n{xi) where is the usual Newtonian potential, 
computed with k as Newton constant. From [26] we get — = (1 + 2. 5 ± 70 x 
10~^){(fis — '^E+ ■ ■ ■ where ^ps is the gravitational potential at the spacecraft 
position and ipE is, the gravitational potential on Earth. . . . accounts for 
additional effects not related to the gravitational potential. We can ascribe 
the uncertainty of the experiment to ^2, to get the bound: 

141 < 24 X 10"*^ 

This bound is conservative because the Newton constant itself has a larger 
error m\: G = 6.67428 ± 0.00067 x IO'^^t^ 

' = kgs 

In our description of the evolution of the Universe, the value of K2 is not 
important, so we will keep it arbitrary for the time being. 

Friedman-Robertson- Walker(FRW) metric This is the main section 
of the paper. We discuss the equations of motion for the Universe described 
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by the FRW metric. We use spatial curvature equal to zero to agree with 
cosmological observations. 

In this paper we will deal only with a perfect fluid, since rotational and 
translational invariance implies that the energy-momentum tensor of the Uni- 
verse has this form. The energy momentum tensor for a perfect fluid is |8]: 

T^u = pgf^u + ip + p)U^U,, g^'U^U, = -1 (16) 

Then: 

^-^r = per + \{P + p){U^U,r + W'U^-g'^^) (17) 

In this case, assuming flat three dimensional metric: 

-ds^ = dt^ - R{tf {dr"^ + r'^de'^ + sin^ ^^0^ j 
-ds^ = A{t)dt^ - B{t) {dr^ + r^de^ + sin^ 6'#'} 

Using ffm [T4|) . we can check that these are co-mobile coordinates and the 
proper time interval dr for a co- moving clock is just dt, so t is the time 
measured in the rest frame of a co-moving clock. Equations ([31 [T7I) give: 

-RB- -pRB + -R-^R^B - -pR^A + -RR^A = 

-pB -2B- R-^R^B + 2R~^RB + 2R~^RB+ 

p R^A + R^A + 2RRA + 2RAR = (18) 

Einstein's equations are: 

3 (;^i?)' f d^ \ f d Y 2 

^ - Kp ,2R{ — R] + i — R] =-kR^p 



We use the equation of state p = wp, to get, for 7^ — 1 : 



2 ~ /w—l\ 

R = i?ot3(i+») ^ = Swht^^', 



B = Rlkt\ b = + (19) 

° ' 3w + 3 w + 1 ^ ^ 



I2 is a free parameter. 
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Red Shift To make the usual connection between redshift and the scale 
factor, we consider light waves traveling to r = 0, from r = ri, along the r 
direction with fixed 9, (j). Photons moves on a null geodesic of g: 

= -(1 + 4i)c/t2 + (i?2 ^ ^'^B){dr^ + r^de^ + sin^ ed(j)'^) (20) 

So, 

A typical galaxy will have fixed r\,Q\,(^\. If a second wave crest is emitted 
at f = ti + bt\ from r = ri, it will reach r = at to + ^^o? where 



ri 



Therefore, for 5ti,5to small, which is appropiate for light waves, we have: 





Introduce: 



We get:|^ = _ A crucial point is that, according to equation ( lT4l) . (5t 

measure the change in proper time. That is:^ = where i^o is the light 
frequency detected at r = corresponding to a source emission at frequency 
vi. Or in terms of the redshift parameter z, defined as the fractional increase 
of the wavelength A: 

z = S^-l = (23) 
Rih) Ai 

We see that R replaces the usual scale factor R in the computation of z. 

Luminosity distance Let us consider a mirror of radius b that is receiv- 
ing light from a distant source. The photons that reach the mirror are inside 
a cone of half- angle e with origin at the source. 
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Let us compute e.The light path of rays coming from a far away source at 
Xi is given hj x (p) = pfi + xi, p > is a. parameter and h is the direction of 
the hght ray.The path reaches us at x = for p = = ri. So n = —xi + e. 
Since n, xi have modulus l,£=|£|<<lis precisely the angle between —xi 
and n at the source. The impact parameter is the proper distance of the path 
from the origin, when p = \xi\. The proper distance is determined by the 
3-dimensional metric (fTSj) . That is b = R (to) ri9 = R (to) ri£, i.e. e = ^^^^^^ . 

Then the solid angle of the cone is vre^ = ^-^t-tt, where A = irb'^ is the 
proper area of the mirror. The fraction of all isotropically emitted photons 
that reach the mirror is / = ^^^2^(-^ -^i ■ Each photon carries an energy hvi 
at the source and hv^ at the mirror. Photons emitted at intervals 5ti will 
arrive at intervals 5to. We have ^ = j^^, = Therefore the power 

at the mirror is Pq = L^p^f, where L is the luminosity of the source. The 
apparent luminosity is / = = L ^^^^j^ inrfRitof ' Euclidean space, the 
luminosity decreases with distance d according to / = ^;j^.This permits to 

define the luminosity distance: 6?^ = 
write this in terms of the red shift: 



define the luminosity distance:(ii = \/ 4^7 = -^(^o)^;^^- Using (12T]) we can 



Supernova la data 

The supernova la data gives, m (apparent or effective magnitude) as a 
function of z. This is related to distance di by m = M + 5log{-^^). Here 
M is common to all supernova and m changes with d^ alone. 

We compare 6 gravity to General Relativity (GR) with a cosmological 
constant: 

= ifo'(^m,(l + Z)' + (1 - fi™)), =l-n^ 

Notice that A = for w = in (fT9|) . So, it seems that we cannot fit the 
supernova data. However = is not the only component of the Universe. 
The massless particles that decoupled earlier still remain. It means that the 
true w is between ^ < | , but very close to w = 0. So, we will fit the 
data with w = 0.1,0.01,0.001 and see how sensitive the predictions are to 
the value of w. 

Using data from Essence j28]. we notice that R^ test changes very little 
for the chosen sequence of w^s. Each fit determines the best I2 for a given w. 
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In this way we see that I2 scales hke I2 ^ a being independent of w. As 
an approximation to the hmit w = 0, we get: 

m = R{t)^^ (25) 

A/a — t 

renormahzes the derivative of at t = 0. It is not divergent, because 

for t ^ 0, w ^ ^. a is a free parameter determined by the best fit to the 
data. 

Of course, the complete model must include the contribution of normal 
matter(ti; = 0) plus relativistic matter (zi; = |). But, at later times, the data 
should tend to ([25]). 

Let us fit the data to the simple scaling model (l25l) . 

We get: 

ilm = 0.22 ± 0.03, M = 43.29 ± 0.03 , x^iperpoint) = 1.0328, General 
Relativity 

a = 2.21 ± 0.12, M = 43.45 ± 0.06, x^iperpoint) = 1.0327, Deha Gravity 

^-gravity with non- relativistic (NR) matter alone give a fit to the data as 
good as GR with NR matter plus a cosmological constant. 

According to the fit to data, a Big Rip will happen at t = 2.21049 in 
unities of to (today). It is a similar scenario as in [23] . 

Finally, we want to point out that since for t — 0, we have — |, then 
R{t) = R{t). Therefore the accelerated expansion is slower than (1251) when 
we include both matter and radiation in the model. 

Conclusions and Open Problems Delta Gravity agrees with General 
Relativity when T^^ = 0, imposing same boundary conditions for both tensor 
fields. In particular, the causal structure of delta gravity in vacuum is the 
same as in General Relativity, since in this case the action is proportional 
to the geodesic action in GR. 

We recover the Newtonian approximation. 

In a homogeneous and isotropic universe, we get accelerated expansion 
without a cosmological constant or additional scalar fields. 

The computation of PPN(Postnewtonian) parameters is in progress, but 
we do not expect large departures from general relativity, because the new- 
tonian limit is the right one, as explained in section 6. Moreover the inter- 
estellar space has very small matter densities, so 6- gravity must give general 
relativity values for the PPN parameters(See comments after equation(jl])). 
Additionally, please notice that all g contributions are multiplied by the small 
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parameter Kg of the order of 10~^ or less, so they are much suppressed in the 
solar system. 

Stellar evolution will not be changed from its newtonian description, un- 
less density of matter becomes very large. Even at the densities of white 
dwarfs the Poisson equation for the gravitational potential suffices. (See, for 
instance [S], chapter 11.3). 6- gravity implies it, as it is shown in section 6. 
Higher densities which are present in neutron stars may provide new tests of 
^-gravity, since there we have to use the whole non-linear Einstein equations 
and the corresponding 6- gravity equations. But for the inner regions of 
massive stars, data is very scarce. 

Notice that equation fll9p implies that R = R aX the beginning of the 
Universe, when w = ^, corresponding to ultrarelativistic matter. That is, 
the accelerated expansion started at a later time, which is needed if we want 
to recover the observational data of density perturbations and growth of 
structures in the Universe. An earlier acceleration of the expansion would 
prevent the growth of density perturbations. 

Work is in progress to compute the growth of density perturbations, the 
anisotropies in the CMBR, BAO, WL and the evolution of massive stars. 
The comparison of these calculations with the considerable amount of astro- 
nomical data that will be available in the near future will be a very stringent 
test of the present gravitational model. 

It was noticed in [20] that the Hamiltonian of delta models is not bounded 
from below. Phantoms cosmological models [22] , [23] also have this property. 
Although it is not clear whether this problem will subsist in a diffeomorphism 
invariant model as delta gravity or not, we want to mention some ways out 
of the difficulty. 

a) Delta gravity is a gauge theory. Moreover it is diffeomorphism invari- 
ant. Thus the canonical Hamiltonian vanishes identically. It may be possible 
to truncate the Hilbert space, using the BRST formalism, to define a model 
with a Hamiltonian bounded from below. This is a difficult task that goes 
far beyond the present paper, but should be pursued in a future work. 

b) In a supersymmetric model we have H = Q^, where H is the Hamilto- 
nian and Q is the hermitian supersymmetry charge. Thus the Hamiltonian 
is bounded from below. So, we expect that a delta supergravity model has a 
Hamiltonian bounded from below. 
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